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Abstract 

We study the spin-Hall effect in systems with weak cubic Rashba spin-orbit interaction. To this 
end we derive particle and spin diffusion equations and explicit expressions for the spin-current 
tensor in the diffusive regime. We discuss the impact of electric fields on the Green's functions and 
the diffusion equations and establish a relationship between the spin-current tensor and the spin 
diffusion equations for the model under scrutiny. We use our results to calculate the edge spin- 
accumulation in a spin-Hall experiment and show that there is also a new Hall-like contribution to 
the electric current in such systems. 
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I. INTRODUCTION 



At present there is much interest in investigations of spin-charge couphng effects in non- 
magnetic semiconductor heterostructures with intrinsic spin-orbit interaction. Particular 
attention is paid to effects, which permit the electric generation of non-equilibrium magne- 
tizations in such systems, like the spin-Hall effect (see, e.g. Ref. uJ - llll ) or the spin accu- 
mulation induced by an external electric field (see, e.g., Ref. [l^- [18[). In the spin -Hall 
effect a non-homogeneous magnetization is created at the sample boundaries in the presence 
of electric fields. The spin accumulation, by contrast, is a homogeneous non-equilibrium 
magnetization, which is induced by a current in the bulk. 

While the mechanisms leading to spin accumulations in heterostructures with intrinsic 
spin-orbit interaction have been understood for decades the spin-Hall effect is still a rather 
controversial issue. The first papers on spin-Hall effects in heterostructures with intrinsic 
spin-orbit interaction have focused on investigations of spin-current tensors for electron gases 
with linear Rashba spin-orbit interaction^ and for hole gases described by the Luttinger 
model-. The spin-current tensors were defined as the expection value of the operator (-Oj- 
velocity operator, Sk -spin operator, {..,..} -anticommutator) 

jik = ^{vi,Sk}. (1) 

The investigations have shown that an electric field leads to distortions of the spectrum, 
which render the spin-Hall current of clean systems finite, even if their Hamiltonian does not 
explicitly contain Sz (S'^-spin-operator for particles with spin in the perpendicular direction). 
This fact has led the authors of the papers in question to the conclusion that there is a non- 
vanishing spin-Hall effect of considerable magnitude in such systems. However, none of 
them has been able to establish a connection between the spin-current tensor and the non- 
equilibrium magnetization induced at the boundary of a sample. The silent assumption was 
the spin current tensor enters a diffusion like equation, as discussed, e.g., in Ref. j7|. 

Further investigations, however, have shown that the results on the spin-current tensor 
are not robust. In the linear Rashba model the spin-Hall conductivity vanishes already in 
the presence of weak disorde r^i^^i^^ . In p-doped systems, described by Luttinger models, 
the situation seems to be more favorable. In such systems the spin-Hall conductivity is also 
finite in the presence of disorderi^*^. In spite of this fact most investigations have focused 
on systems with linear Rashba spin-orbit interaction, which are technically less demanding. 

The absence of spin-Hall currents in the linear Rashba model in the presence of weak 
disorder, however, has not remained the only problem casting doubts on the predictions of 
spin-Hall effect. More pressing is the question of the relevance of the approach itself^^^. It 
has already been shown in Ref.^J that the spin-current tensor in the Rashba model is finite, 
even in equilibrium. Therefore, it does not characterize macroscopic transport phenomena. 
Diffusion equations, which take into account also spin-charge coupling effects, have been 
derived in a number of papers for systems with weak Rashba interactio n^^i^^i^^i^'''i^^i^^i^°i^'^ . 
However, they do not only contain the divergence of the spin-current tensor but also addi- 
tional derivative terms. The investigation of these terms leads to the conclusion that the 
spin-current tensor in the linear Rashba model is not a quantity, which is continuous at 
boundaries^!, as stressed also in Ref. 2^. In fact, the solution to the diffusion equations. 



which is derived under the assumption that the spin-current tensor vanishes at hard-wall 
boundaries, is not physical. It violates time reversal symmetry- Responsible for this fact 
are just those terms, which are also present in equilibrium. Alternative definitions for the 
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spin-current tensor have been suggested in a number of papers (see, e.g., Refs. |32ll33j ) 
to circumvent these problems. However, although the authors have focused on observable 
quantities they have failed to establish a connection between the non-equilibrium magneti- 
zation at the boundary of a sample and their definition. Moreover, their predictions differ 
rather strongly from those, which have been derived on the basis of Eq.([Q). The existence 
of the spin-Hall effect in clean samples has meanwhile been established experimentally^ in 
spite of these problems. The predictions of Ref.yj have essentially proven to be correct. The 
question, however, how to describe this effect theoretically, is still remaining. 

Recently, some investigations have been published, which address the question how uni- 
versal the results to the linear Rashba model are^^iSi^i. The results of these investigations 
show that already deviations from the quadratic dispersion relation of the unperturbed band 
spectrum'^'* or deviations from the linear dependence of the Rashba spin-orbit interaction^^'^— 
lead to a non-vanishing spin-Hall current in the bulk even in the presence of disorder. To 
obtain the latter result the authors of the Refs. l3] and [3^ use a cubic Rashba model. The 
cubic Rashba model, which describes the heavy holes in III-V semiconductor quantum wells, 
has been introduced in the Refs. j36il37l | . It differs from the linear Rashba model in that 
the interaction depends cubically on the momentum. Thus, the cubic Rashba model offers 
the opportunity to study spin-Hall transport in semiconductor heterostructures with intrin- 
sic spin-orbit interaction in a framework, which keeps the simplicity of the original Rashba 
model. Here we use this opportunity to address the question, how the non-vanishing spin- 
Hall current manifests themselves in the spin diffusion equations, how these equations differ 
from those in the linear Rashba model, and whether the spin-current tensor is the ap- 
propriate quantity, which has to be considered, or whether alternative definitions for the 
spin-current tensor are more appropriate, as claimed, e.g., in the Refs.js^js^. To answer 
these questions we derive particle diffusion equations, spin diffusion equations and explicit 
expressions for the spin-current tensor, discuss the impact of external fields on the diffusion 
process and use our results to calculate the magnitude of the magnetization for systems with 
weak Rashba spin-orbit interaction. Doing so, we establish an explicit relationship between 
the spin-current tensor (^, the diffusion equations, and the non-equilibrium magnetization 
at the boundary in the presence of an electric field. Thus our results show that the spin- 
current tensor (Q) determines indeed the physics in the cubic Rashba model. We would like 
to mention that a paper has recently been published, in which the spin-Hall effect has been 
investigated numerically in the cubic Rashba modeR The authors of Ref . [9] focus on sys- 
tems with strong Rashba interaction, in which the energy level splitting at the Fermi surface 
is large compared to the disorder energy. In contrast to them we present analytical calcu- 
lations and focus on systems in the opposite regime. While the numerical results of Ref. 
apply to systems with strong spin-orbit interaction, in which the energy level splitting at 
the Fermi energy due to the spin-orbit interaction is large compared to the disorder energy, 
we focus on systems with weak spin-orbit interaction, in which the energy level splitting at 
the Fermi energy is small compared to the disorder energy. 



II. BASIC EQUATIONS 

We consider the holes in the lowest heavy-hole subband of a two-dimensional hole gas with 
Rashba spin-orbit interaction. Their motion is described by the Hamiltonian^j^^i^^iSL^Si^ 

H=^ + (P' ^+ - Vl^-) + V{r). (2) 
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Here p is the momentum operator, m is the effective mass, p± = px ± ipy, a± = ax i: iCy 
and is a constant characterizing the strength of the spin-orbit interaction. The potential 
V{r) is a random white noise potential with Gaussian statistics and correlation function 

{Vir)V{r')) = -^6{r-r'), (3) 

where u = m/27i?i^ is the density of states per spin and r is the scattering time. The vectors 
Bx and By are the unit vectors in the two-dimensional plane, is a unit vector perpendicular 
to the plane. The indices 1,2 and 3 replace the indices x,y and z in sums. Double indices 
have to be summed over in all formulas. Sums over Latin indices run from 1-3, sums over 
Greek indices from 0-3. 

The Hamiltonian (j2]) describes particles with single particle energies E±{k) = ± A^, 
where = h'^k'^/2m, hk is the momentum of the particles, and = \N\h^k^ is the energy 
level splitting due to the Rashba spin-orbit interaction. The spectrum is unstable at large k, 
since A^. oc k^. This instability is an artifact, which results from the fact that the effective 
Hamilton (0) is only derived for systems with small k by means of perturbation theory. 
Therefore, the Hamiltonian Q can only be used for low hole densities. From the practical 
point of view these are also those typical achieved in experiments^^,. The investigation of 
transport processes on the basis of the Hamiltonian however, requires also the calculation 
of integrals with infinite upper bound. To render them finite we assume that A^ depends 
only cubically on k in the region of interest, but increases slower than k^ for k —>■ oo. In this 
case we can use the same methods for the investigation of the system, which have also been 
used for systems with linear Rashba spin-orbit interaction (see, e.g., Refs. [24,..25..26i27ll4^ ) . 
In particular we can use the self-consistent Born approximation for the description of the 
single-particle excitations. The retarded Green's function takes the form 

G^ik\E) = gl{k\E)g^{k\E){g],\k\E)a, + i^{pla^ - pla^)) (4) 
in this approximation, where 

(To = 1 is the unit matrix and gR[k\E) = g^{k\E)\^=Q. The advanced Green's function is 
obtained from Eq.(@]) by hermitian conjugation. 



III. SPIN DIFFUSION IN THE ABSENCE OF EXTERNAL FIELDS 



A. Particle and spin diffusion equations 



To derive particle and spin diffusion equations we use the ladder approximation. The 
equation of motion for the Fourier-Laplace transform of the particle density n{K,\E,s) = 
Sq{k\E,s) and the Fourier-Laplace transform of the spin density S{k,\E,s) = Si{K,\E, s)ei 
takes the formSiSLM 

n r n if K K 

K7 7^trKG^(fc + 2 1^ + ^ns)a,G^ik - -\E))]S,iK\E, s) = So,i>^\E) (6) 
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in this case. Here s is the frequency describing the evolution (it corresponds to a Laplace 
transformation with respect to time), tr symbolizes the trace, and Sorjif^lE) is the initial 
condition. To simplify Eq. © we use the hydrodynamic expansion, viz we expand the kernel 
of Eq.® with respect to k,, N and s and take into account only contributions, which are 
at most of second order in k and of first order in s. In this case we obtain the following 
diffusion equations for systems with weak Rashba spin-orbit interaction after inverse Laplace 
and Fourier transformation: 

dtn - DAn = 0, (7) 

dtS + n-S- DAS = 0. (8) 

Here D = Er/m and f2 is a symmetrical tensor of second rank with matrix elements 
Oifc = VL5ik{l + 5a), where 

Vt = 16 ^E^T. 9) 

In deriving these equations we have ignored terms of the order N'^'VfSi, which yield cor- 
rections to the diffusion coefficient. These corrections are small with respect to the ratio 
(Afc|fj.=e^/(^/2r))^ and thus negligible in systems with weak Rashba spin-orbit interaction. 

It is instructive to compare these equations with those, which are obtained for systems 
with linear Rashba spin-orbit interaction in the same approximation. If we replace the cubic 
spin-orbit interaction in Eq.((2} by a linear interaction, according to 

(P' - pV-) ^{N,cTxp), (10) 
where N = Ne^, we obtain the equations^iiS^iSLSLS^ 

dtn - DAn - r]r(V, N x S) = 0. (11) 

dtS + n-iS- Sa) - DAS - u},{N X V) X (5 - Sa) = 0. (12) 

Here 

Sa = -tN X Vn (13) 

is the non-equilibrium magnetization induced by a density gradient and Ug is a transport 
coefficient proportional to D (Note that, although we use the character to characterize the 
strength of the Rashba spin-orbit interaction both for systems with linear and cubic Rashba 
interaction this quantity is different in both systems and possesses different dimension. The 
structure of the tensor O in the linear model agrees with that in the cubic model. The 
dependence of the quantity Q on the microscopic details, however, is also different in both 
models.). 

Comparing the Eqs.© and (jH)) with the Eqs. pi|) and p2|) we note the following. Whereas 
the spin-orbit coupling gives rise to a coupling between spin and charge in the linear model 
there is no coupling between spin and charge in the diffusion equations for cubic systems 
at all. There is neither a spin accumulation {Sa = 0) nor a spin-galvanic current. This 
fact is very surprising, since we would have expected that the structure of the diffusion 
equations is little affected by a simple change in the power of the momentum in the Rashba 
spin-orbit interaction. We would like to mention that we have also checked this conclusion 
by calculating the spin accumulation in the Kubo-formalism. The results of this calculation 
are in line with those obtained from the diffusion equation and therefore not presented 
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here. They also show that the dc-spin accumulation is zero in the linear approximation with 
respect to the electric field for systems with weak cubic Rashba spin-orbit interaction. 

The diffusion equations (jH)) and ()12|) . however, differ not only in the spin-charge coupling 
terms. Comparing both equations further we also note that there is no term linear in the 
derivatives in the spin-diffusion equation (jl2j) . neither in the linear approximation with 
respect to N nor in higher order. This fact has far reaching consequences. The coefficient 
Us in Eq. (|12j) leads to a precession of the magnetization around an axis, which is parallel to 
N X F,i{ an electric field F is switched on. This precession can be used for the manipulation 
of the electron spin (see, e.g., Refs. |27ll30ll40lj4l|j42l ]). The corresponding effective Zeeman 
field is linear in F. In the cubic Rashba-model by contrast, there is no effective Zeeman-field 
in linear approximation with respect to F. 

The huge difference between the structure of the diffusion equations raises the question 
why the physics of both systems is drastically different. There is no answer to the question 
yet. On the technical level, however, the difference results from the angle averaging. The 
cubic Rashba-interaction is proportional to p\ = {hk±)^. The calculation of terms of first 
order in the derivatives, e.g. of spin-charge coupling terms, reduces to integrals of the type 

which are zero for arbitrary functions f{k). 



B. The spin-current tensor 

The surprising difference between the structure of the spin diffusion equations for systems 
with linear Rashba spin-orbit interaction and that for systems with cubic Rashba spin-orbit 
interaction raises the question, where spin-charge coupling effects manifest themselves. To 
investigate this question further we investigate the spin-current tensor. The spin-current 
tensor is defined as the expection value Jik{r\E, s) of the operator 

Mr) = l{jik,Sir-r)} (14) 

where ^ 

3ik = ^Wi,Vk}, (15) 

f is the position operator and Vk is the velocity operator. To calculate the expection value 
we restrict the consideration to time scales, which are large compared to the momentum 
relaxation time. In this case we can use the following relationship between the spin-current 
tensor and the Green's function o^^i^^ : 

Mr\E, s) = f dr,drMG^{ri - r2|^)4(r)G'^(r2 - n\E)a^)S^{n\E, s). (16) 

Here the Green's functions are those in the position representation. To simplify the integrals 
we use again the hydrodynamic expansion and obtain the equations 

Ji, = -D\/^S,, J,y = -DViSy (17) 
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and 

3 h 

Ji, = -DViS, + --TVt{e, X V)in. (18) 
2m 

Eq. (fT^ generalizes the existing results^ on the spin-Hall current to non-homogeneous sys- 
tems. To compare Eq. (fTHj) with the result of Ref . [35| for the spin-Hall current of a homoge- 
neous system in the presence of an electric field we use the fact that the impact of the field 
on the spin-Hall current can be mimed by means of a concentration gradient in the ohmic 
approximation^^. To mime the field we use a particle density of the form 

n{r\E,s) =no-2v{F,r) (19) 

where F = —eS (e is the electric charge and S is the electric field) and uq is a constant 
density. Doing so, we obtain 

Ji, = -3n-Ti^{e, X F)i. (20) 

7Yt 



This result coincides with that of Ref. j35|] . It shows that there is a non- vanishing spin-Hall 
current in this system. The crucial question is, whether this fact gives also rise to a non- 
vanishing spin-Hall effect, viz to a non-equilibrium spin accumulation at the boundaries of 
the sample. 

The equations ()17p and (jlSp differ strongly from those obtained for systems with lin- 
ear Rashba-interaction^L^i. The spin-current tensor for such systems is even non-zero in 



equilibrium (see, e.g. Refs. |22ll23| ). This fact manifests itself in terms containing the par- 
ticle density without derivative a^^ . Such terms are absent in systems with weak cubic 
Rashba-interaction. Therefore, the spin-current tensor describes deviations from equilibrium 
in systems with weak cubic Rashba-interaction and thus characterizes transport phenomena. 
In contrast to the spin diffusion equation, however, the spin-current tensor couples also spin 
with charge. The coupling is achieved by the second term on the right hand side (rhs) of 
Eq. lfTHj) . Since such spin-charge coupling terms are absent in the diffusion equation there is 
the question, whether there is really a connection between the spin diffusion equation and 
the spin-current tensor. Such a connection has not been established for systems with linear 
Rashba-interaction yetS^i. 

From the phenomenological point of view we expect that the spin diffusion equation can 
be written in the forn>^ 

dtSk + ClkiSi + Vj,k = 0. (21) 

Here Jik is a tensor of currents, which we call transport currents. The transport currents 
determine the boundary conditions needed to find solutions to the diffusion equation. The 
physical solutions are those with the property that the transport currents are continuous at 
the boundary in the presence of a boundary^^. To find the transport currents there is no 
other way than to read off all terms, which are under the sign of the derivative. However, 
since Eq. ()21|) contains only the divergence of the spin-current tensor it is impossible to 
retrieve terms proportional to a curl unambiguously. 

The spin-charge coupling term in Eq. lfTHjl has just just the structure of a curl. Due to 
this fact the relationship 

V,J,fc = -DASk (22) 
holds, so Eq.® can also be written in the form 

dtSk + CtkiSi + V^J^k = 0. (23) 
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This fact suggests that the spin-current tensor coincides with the tensor of transport currents. 
Spin-charge couphng effects would enter the diffusion equation via the boundary condition 
if this assumption were correct. 

Close to equilibrium, however, there is no way to assure that the spin-current tensor really 
coincides with the tensor of transport currents, since the terms proportional to a curl can 
not be retrieved from the diffusion equation, as mentioned before. Thus, there is no reason 
to believe that the spin-current tensor really determines the boundary conditions. There is 
no generally accepted method to avoid this problem. However, in a recent papep^I we have 
noted the following fact: we can get further insight into the structure of the terms in question 
by switching on an additional external electric field F. The transport coefficients become 
position dependent in this case, since they depend on the kinetic energy fir = E — {F, r) of 
the particles, which is different for particles with fixed total energy E at different positions 
r. Due to this fact the derivatives in the diffusion equation have also to act on the transport 
coefficients and thus, they produce additional terms. If the assumption that the spin-current 
tensor coincides with the tensor of transport currents is correct, these terms are just those, 
which result from derivatives of the spin-charge coupling term on the rhs of Eq. (jl8p . Thus, 
this observation is helping us to identify unambiguously all those terms in the tensor of 
transport currents, which have the structure 

uj{fir){ez X V)^S^{r\E,s), 

where uj{fir) is an arbitrary transport coefficient. Only terms of the form 

(e^ X V)i{uj{nr)S^{r\E,s)) 

can not be retrieved. 



IV. SPIN DIFFUSION IN THE PRESENCE OF AN ELECTRIC FIELD 
A. Diffusion in the presence of an electric field 

To generalize the diffusion equations to the presence of an electric field we use the sym- 
metry properties of the Green's functions. The Green's functions satisfy the property 

G{r + a,r' + a\E) = G{r, r'\E - (F, a)) (24) 

in the presence of a constant field F, where G{r,r'\E) is either the retarded or the advanced 
Green's function in the position representation and a is an arbitrary vector. This fact permits 
the introduction of a new function g, according to the relationshipi^^ 

^^'^' "^''^^ = / J^2^'^'''^"''^~9(f'M, (25) 

where R = [r + The knowledge of the function g is sufficient to derive the diffusion 

equation. In the presence of the field it takes the forn>^ 

M^^r., s) + V^} + \v.^^,{^^r)Vi]Ss{r\E, s) = So,{r\E). (26) 
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where {...,...} is the anticommutator. The coefficients a, b and c in this equation are 
related to the function g by the relationships 

a^5(/ir, s) = 5^5 - J j^^tT[a^g^{k\iJ,r + ihs)asg^{k\ij,r)], (27) 
^■ysM = J -^tT[a^g''ik\fir)cT5tVjg^ik\fir)] (28) 



and 



AnvT J (2-K 
Ti f dk 

= ^ J j^tT[a^g''{k\fir.)crsViVkg^{k\fir)]. (29) 

The function g can be calculated perturbatively. It depends only on in the absence of 
the spin-orbit interaction and reduces to the equilibrium Green's function in such systems 
in the ohmic approximation. The coefficients a, b and c are the same in this case as those in 
the absence of the field. Therefore, the impact of electric fields can be reduced to position 
dependent shifts of the reference point of the kinetic energy, which can be taken into account 
by a simple substitution of the type V ^ V — Fd^ in the equilibrium diffusion equation, 
where e is the kinetic energy. 

However, in the presence of the spin-orbit interaction this concept proves to be inadequate. 
The impact of the field on transport processes can not be reduced to simple shifts of the 
reference point of the kinetic energy in such systems, as claimed, e.g., in the Refs.j^ and 
(28| . In such systems the field leads both to a position dependent shift of the reference 
point of the kinetic energy and to distortions of the particle spectrum. These distortions 
are responsible for the spin-Hall effect, as discussed in Ref.0]. To take into account the 
distortions we restrict the consideration to the ohmic approximation. In this case we can 
write the function g in the form 

~gik\E) = G{k\E) + ~g^'\k\E), (30) 

where g^^\k\E) is the linear contribution to g with respect to F and the function G{k\E) 
is the equilibrium Green's function. The function g^^^ can be calculated perturbatively. For 
the linear Rashba-model we obtain g^^^ oc N'^{F, crxk) (see Eq.(lO) of Ref. ji^), as expected 
of the spin-Hall term. The cubic Rashba-model yields 

g^'\k\E) = t-a^N\g^{k\E)g4k\E)y{plplF^-plplF^), (31) 

where F± = F^ i: iFy. If we take into account this term in calculating the coefficients in 
Eq. p6|) we obtain the equation 

dtSk{r\E, t) + V^ki{^lr)S^{r\E, t) + Vi J.fc(r|E, t) = (32) 

in the ohmic approximation, where ri(/ir) = ri|_B=^^ and the spin-current tensor is given by 
the Eqs.dni) and ^ with Q. fi(/i^) and D D{fir) = D\E=f^,. In Eq.^ we have also 
included the variables the densities depend on to stress that the densities are functions of 
the total energy E, of the particle position r and of t in the presence of a field. Below we 
omit the variables again to avoid cluttering the notation. 

At this point it is important to note the following: whereas Eq.(jHJ can only be written 
in the form due to the identity (|^. Eq.(j221) contains also derivatives of spin-charge 
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coupling coefficients, since the equality (j221) does not hold in the presence of electric fields. 
The spin-charge coupling term (second term on the rhs of Eq. ljlSp ) arises naturally in the 
diffusion equation. It does not reduce to the divergence of a curl. Consequently, we can 
truly claim that the spin-charge coupling terms enter the diffusion equation and thus, that 
the spin-current tensor coincides with the tensor of transport currents. The second point 
we would like to draw attention to is the following: the spin-charge coupling terms would 
be absent if we had ignored the distortion of the particle spectrum due to the electric field. 

Eq.(j221) couples the spin density with the charge density. The diffusion equation for the 
charge density takes the form 

dtn{r\E, t) + divj{r\E, t), (33) 

where 

j{r\E, t) = -D{fir)Vn{r\E, t) - — -n(/x,)e, x VS,{r\E, t). (34) 

2 m 

Eq.()34|) sets the situation again apart from that in the linear Rashba-model. Homogeneous 
non-equilibrium magnetizations in the plane of the 2-d electron gas lead to electric currents 
in the linear Rashba-model in the absence of external fields (the spin-galvanic currents). 
In the cubic Rashba model, by contrast, non-equilibrium magnetizations can only induce 
electric currents in the absence of external fields if they are inhomogeneous and perpendicular 
to the 2-d plane of the 2-d electron gas. A homogeneous non-equilibrium magnetization 
perpendicular to the 2-d plane can only induce a current in the presence of a field, as 
discussed further below. 

The Eqs. ()32|) - (jM|) do not have the structure of conventional diffusion equations yet. The 
reason for this fact is that these equations are written down for densities, which depend on 
the total energy of the particles. To obtain conventional diffusion equations we distinguish 
between the position dependence of the densities, which characterizes the shape of the packet, 
and the position dependence of the reference point of the kinetic energy, viz we write the 
densities in the form n{r\E,t) n(r|/i^,t), S{r\E,t) — > S{r\fj,r,t) and replace fir by e. 
The position dependence of the density n{r\e,t) describes the shape of a particle packet of 
particles with the same kinetic energy e, which is measured from the tilted bottom of the 
band. Likewise for the density S{r\e,t). Doing so, we find that the derivatives transform 
according to the rule V — V — Fd^. Using this method, we obtain 

dtSk + nk^ie)S^ + WiJik = -d,j,i, (35) 

where j^^ is the energy relaxation current for spins, 

J,, = -D{e)VA + D{e)FAS., J^y = -DV^Sy + D{e)FASy, (36) 

and 

J,, = -D{e)V,S-, + D{e)FAS. + --rf^(e)(e, x V),n - --rf^(e)(e, x F)An. (37) 

2 m 2 m 

The second term on the rhs of the Eqs. ()36|1 and p7jl might look unfamiliar. However, in 
looking on these equations one has to take into account that these equations are spectral 
diffusion equations. To obtain diffusion equations, say for particles at the Fermi surface, 
one has to integrate these equations with respect to e. In this case the first term on the rhs 
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of Eq.(jnni) vanishes and the second term on the rhs of the Eqs. fjH^ and (|H7j) yields just the 
mobihty times the field, as expected. 

The electric current takes the following form 

3^ r 3S r 

j = -D{t)Vn + D{t)Fd,n n(e)e, x V^, + Vt{t)e, x Fd,S,. (38) 

2 m 2 m 

and the diffusion equation (fH^ yields 

dtn + divj = -(9eje, (39) 

where is the energy relaxation current for particles. Note that, the fourth term on the 
rhs of Eq. ()38|l yields also an electric current in an electric field if there is a homogeneous 
magnetization in z-direction. This current is a Hall current, i.e., it is transverse to F. 



B. Spin-Hall effect 

To investigate the spin-Hall effect we focus on a system in the half-plane y > 0. We 
assume that the system is translation invariant in x-direction [d/dx = 0), that the electric 
field is applied in x-direction and that there is an additional scattering process, which keeps 
the system in energy equilibrium. Thus, n{r\e,t) = no6{eF — e), where ep is the Fermi- 
energy and no = 2i/. A glance on the spin diffusion equations (jH^j) - (jHTj) reveals that a 
non-equilibrium magnetization at the Fermi surface is created in this case. Therefore, we 
can solve the spin diffusion equation with the ansatz S{r\e,t) = S{r\t)6{e — ep). To obtain 
an equation for S{r\t) we restrict the consideration to the stationary limit and integrate the 
spin diffusion equations with respect to e. Doing so, we obtain the equation 

2n{eF)S,-D{tF)^S, = Q, (40) 
which has to be solved subjected to the boundary conditions 

= Jy,\y=o = -D{eF) — S, + —n{tF)TF,no (41) 
and Sz = finite for ?/ — oo. The unique solution is 

3 



2 \ 2D{eF)m' 



TnoF,exp{-J2n{eF)/D{eF)y). (42) 



Thus, the spin-Hall current in the bulk leads to a non-equilibrium magnetization at the 
boundary, as it is also the case in systems with strong Rashba spin-orbit interaction^. Since 
the exponent depends only on the ratio Q{eF)/ D^ep) the penetration depth of the spin 
accumulation is independent of disorder in systems with weak spin-orbit interaction. It 
decreases with increasing particle number. The magnitude of the edge spin-accumulation 
increases with increasing ep since ^2(6^) oc . This sets the situation apart from that 
in nearly clean systems. There both the magnitude of the edge spin-accumulation and the 
penetration depth decrease with increasing particle number—. 

The disorder enters only in the preexponential factor via the relaxation time. This fac- 
tor leads to a reduction of the magnitude of the magnetization with increasing amount of 
disorder. 
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C. Hall-like current 



Spin-charge coupling effects manifest themselves also in the electric current in optical ex- 
periments. To investigate the magnitude of the Hall-like current, we assume that a constant 
non-equilibrium magnetization of the type S{r\e,t) = S^e^S^e — ep) is created optically by 
means of optical orientation. Such a density leads to a Hall-like current of the form 

Sj = --—e, X FS, ' . (43) 
2 m aep 

Surprisingly this current depends only quadratically on A^, as it is also the case for the spin- 
Hall effect. It shares this property with the Hall current, which has recently been discussed 
in Ref . [46^ . We speculate that the physics is similar. 



V. CONCLUSIONS 

We have studied the diffusion equations for heavy holes in semiconductor quantum wells 
with Rashba spin-orbit interaction. The holes are described by a cubic Rashba-model. 
In our investigation we have restricted the consideration to systems with weak Rashba 
interaction. These systems have the property that the energy level splitting A due to the 
Rashba interaction at the Fermi energy is small compared to the disorder energy h/2T. 

At first glance the cubic Rashba model differs from the conventional linear Rashba model 
only in that the interaction does not depend linearly but cubically on the momentum. This 
small difference, however, leads to drastic changes in the structure of the spin diffusion equa- 
tions. The particle and spin diffusion equations reduce to conventional diffusion equations 
in the absence of external fields. There are no spin-charge coupling effects in the diffusion 
equations themselves in the absence of external fields. There is neither a spin-accumulation, 
which can be induced by concentration gradient, nor a spin-galvanic current in homoge- 
neous non-equilibrium systems. This sets the situation strongly apart from that for systems 
with linear Rashba-interaction, in which spin-charge coupling terms occur in the particle 
and spin diffusion equations, which are related to the existence of these phenomena. If we 
switch on an external field the situation changes. Spin-charge coupling terms occur also 
explicitly in the diffusion equation, which lead to a non-equilibrium current in the presence 
of a perpendicular non-equilibrium magnetization and to spin-Hall effect. Our investigation 
shows that both effects depend only on the magnitude of the Rashba interaction strength. 
Explicit formulas for the magnitude of the spin-accumulation at the boundary and for the 
magnitude of the non-equilibrium current under steady state illumination are given in the 
text. 

The investigation of the structure of the spin-diffusion equations has also yield further in- 
formations on the physics of spin-diffusion processes in such systems. It has shown that the 
spin-current tensor coincides with the tensor of transport currents in the diffusion equation 
(see Eq. (|2ip for the definition of tensor of transport currents). Therefore, the spin-current 
tensor determines the boundary conditions needed to find solutions to the spin diffusion 
equations. Due to this fact spin-charge coupling effects enter the solution to the diffusion 
equations even in the absence of external fields, although they are absent in the diffusion 
equations themselves in the absence of fields. This fact is remarkable and sets the situa- 
tion again strongly apart from that in the linear model. In the linear model there is no 
relationship between the tensor of transport currents and the spin-current tensor—. 
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The huge differences between the results for the hnear and the cubic model raise the 
question, why the physics in the cubic Rashba model differs so strongly from that in the 
linear model. We have no clear answer to the question yet but speculate that the difference 
is caused by the fact that the linear spin-orbit coupling has nearly the structure of a pure 
gauge. In fact, to take into account the spin-orbit interaction in the linear model we only 
have to replace the derivatives by covariant derivatives in the diffusion equations, as long as 
spin-charge coupling effects are ignored^. 
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